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Coherent control of a system involves steering an interaction to a final coherent state by con¬ 
trolling the phase of an applied field. Plasmas support coherent wave structures that can be 
generated by intense laser fields. Here, we demonstrate the coherent control of plasma dy¬ 
namics in a laser wakefield electron acceleration experiment. A genetic algorithm is imple¬ 
mented using a deformable mirror with the electron beam signal as feedback, which allows a 
heuristic search for the optimal wavefront under laser-plasma conditions that is not known a 
priori. We are able to improve both the electron beam charge and angular distribution by an 
order of magnitude. These improvements do not simply correlate with having the ‘best’ focal 
spot, since the highest quality vacuum focal spot produces a greatly inferior electron beam, 
but instead correspond to the particular laser phase front that steers the plasma wave to a 
final state with optimal accelerating fields. 

1 Introduction 

The concept of coherent control — precise measurement or determination of a process through 
control of the phase of an applied oscillating field — has been applied to many different systems, 
including quantum dynamics', trapped atomic ions^, chemical reactions^. Cooper pairs"', quantum 
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dots^’® and THz generation^ to name but a few. A plasma wave is a coherent and deterministically 
evolving structure that can be generated by the interaction of laser light with plasma. It is therefore 
natural to assume that coherent control techniques may also be applied to plasma waves. Plasma 
waves produced by high power lasers have been studied intensively for their numerous applica¬ 
tions, such as the production of ultrashort pulses by plasma wave compression^’ generation of 
extremely high power pulses by Raman amplification^^’for inertial confinement fusion ignition 
schemes^^’as well as for fundamental scientific investigations. In particular, laser wakefield 
acceleration of ultra-relativistic electron beams'^“^^, has been a successful method for accelerat¬ 
ing electrons to relativistic energies over a very short distance. In laser wakefield acceleration, 
an electron bunch ‘surfs’ on the electron plasma wave generated by an intense laser and gains a 
large amount of energy. The accelerating electric field strength that the plasma wave can support 
can be many orders of magnitude higher than that of a conventional accelerator, which makes 
laser wakefield acceleration an exciting prospect as an advanced accelerator concept. However, 
although highly competitive in terms of accelerating gradient, beams from laser wakefield accel¬ 
erator experiments are currently inferior to conventional accelerators in terms of other important 
characteristics, such as energy spread and stability. In addition, due to constraints in laser wake¬ 
field technology, experimental demonstrations have predominantly been performed in single shot 
operation, far below the kHz-MHz repetition rates of conventional accelerators. 

In recent years, deformable mirror adaptive optical systems have been successfully imple¬ 
mented in high intensity laser experiments to increase the peak laser intensity by improving the 
beam focusability, especially in systems using high numerical aperture optics. The shape of the 
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deformable mirror is generally determined in a closed loop where either a direct measurement of 
the wavefront is performed^^ or some nonlinear optical signal^®’is used as feedback in an iter¬ 
ative algorithm. The objective of adaptive optics has largely been optimization of the laser focal 
shape to a near diffraction-limited spot, thus producing the highest possible intensity. Adaptive 
optics can also be useful for certain focal profile shaping^^’^^, optimization of a laser machining 
process^"^ or harmonic generation^^’^®. 

In the following, we demonstrate that orders of magnitude improvement to electron beam 
properties from a laser wakefield accelerator operating at kHz repetition rate can be made, through 
the use of a genetic algorithm coupled to a deformable mirror adaptive optical system to coherently 
control the plasma wave formation. The electron image from a scintillator screen was processed 
and used in the fitness function as feedback for the genetic algorithm. Using this method, we 
were able to improve the beam properties significantly. This result was not simply due to an 
improvement in focal quality since a laser pulse with the ‘best’ (highest intensity/lowest M^) focus 
in vacuum produced a greatly inferior electron beam compared with a laser pulse optimized using 
the electron beam properties themselves. It was found that the focal spot optimized for electron 
beam production had pronounced intensity ‘wings’. Modifications to the phase front of the tightly 
focusing laser alter the light propagation, which experiences strong optical nonlinearities in the 
plasma, and therefore affect the plasma wave dynamics in a complex but deterministic manner. 
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2 Results 


Experimental setup and procedure The experiment was performed using the relativistic Lambda- 
Cubed (A^) laser system (see Methods). The output laser beam was reflected from a deformable 
mirror and focused onto a free-flowing argon gas plume to produce an electron beam by laser 
Wakefield acceleration (see Methods) at 500 Hz. Electrons were measured using a scintillating 
screen imaged onto a lens-coupled CCD camera. The experimental setup is shown schematically 
in Fig. 1. 

We first implemented a genetic algorithm for laser focus optimization using the second- 
harmonic signal generated from a beta barium borate (/5-BBO) crystal (setup A in Fig. 1). The 
laser spot was optimized such that highest peak intensity is achieved when the second harmonic 
generation is strongest. Subsequently, we modified the fitness function to use a figure of merit 
(FOM, refer to equation 1 in Methods) from the electron scintillation data, calculating the inverse 
distance weighting (with power parameter n) to a single point fq for all pixel intensities within an 
electron image. The pixel of the optimization point Fq was dynamically adjusted during the genetic 
algorithm to concentrate all electron signal to the peak location of the charge distribution during 
each generation. The genetic algorithm was initialized using a ‘flat’ mirror shape with 30 V for all 
actuators to allow immediate deformation in both directions. 

Optimization of the electron spatial profile For comparison, electron beams produced by the 
‘best’ laser focus (by optimizing the intensity) and the initial mirror shape at 30 V are shown in 
Fig. 2a and b respectively. The optimized electron beam profiles are shown in Fig. 2e-j for various 
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weighting parameters, n. The genetic algorithm converged to the best electron beam using n = 8 
in terms of beam divergence and peak charge density. The peak charge density was increased 
by a factor of 20 compared to the initial electron beam profile before optimization (see Fig. 2d). 
The optimized electron profile is highly stable and collimated, with a full-width-at-half-maximum 
(FWHM) divergence of Aa: = 7.4±0.6 mrad and Ay = 12.8±1.4 mrad. The shot-to-shot pointing 
(defined by the centroid position) fluctuation of the electron beam is less than 1 mrad (root mean 
square, r.m.s.). The integrated charge was increased by more than two-fold from the electron beams 
generated by a laser focus of highest intensity. The high repetition rate and real-time diagnostics 
permit implementation of the algorithm within a practical time frame using a standard personal 
computer. Typical optimization takes only a few minutes (~40 iterations) to reach convergence 
(see Fig. 2c). 

The second harmonic optimization^® generates a near-diffraction-limited focal spot as shown 
in Fig. 3a for the far-field laser intensity profile in vacuum. In Figs. 3a and b, we compare the 
transverse intensity distribution within the focal region over the length scale of the gas jet. The 
laser profile (Fig. 3b) that produces the best electron beam exhibits several low intensity side lobes 
around the central peak, and has a peak intensity about half that of the optimized focus. The 
complex laser profiles appear to have a very dramatic effect on the structure of the plasma waves 
produced and consequently the electron beam profile. Fig. 3c shows the relative wavefront change 
recorded by a Shack-Hartmann wavefront sensor. Calculation using the reconstructed wavefront 
gives a Strehl ratio of about 0.5, which is in agreement with the far-field intensity measurement. 
This small wavefront modification of the driver pulse can lead to a significant improvement in 
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the electron beam properties through the relativistic nonlinear optics of the plasma. The relative 
position between the focal plane and the center of the gas flow was controlled by moving the 
nozzle. Scanning the gas nozzle both before and after genetic algorithm confirms the optimal focal 
position does not change, excluding the possibility that the improvement may be due to optimizing 
the focal position. 

Control of energy distribution Furthermore, we extended the genetic algorithm optimization to 
control the electron energy distribution. Through control of the light propagation, the plasma wave 
amplitude will be affected and therefore also the strength of the accelerating gradient. Hence, we 
can expect to be able to modify the energy spectrum. A high resolution energy analyzer using a 
dipole magnet pair was used to obtain the electron energy spectrum as the electrons were dispersed 
in the horizontal plane in the magnetic field. A 150 /um pinhole was placed 2.2 cm from the electron 
source to improve the energy resolution of the spectrometer. The schematic setup is shown in 
Fig. 4a. The energy resolution limited by the entrance pinhole and transverse emittance of the 
beam is estimated to be 2 keV for the energy range of measurement. 

Three rectangular masks are set in the low-, mid- and high-energy region on the dispersed 
data, namely masks I, II, and III in Fig. 4b. We employed a fitness function (see Methods) to 
preferentially maximize the total counts inside the mask. Raw spectra from the genetic algorithm 
optimization are displayed in Fig. 4b, showing that the brightest part has shifted congruently. The 
resulting spectra have mean energies of 89 keV, 95 keV and 98 keV respectively for masks I, II, 
and III, noting that they do not fall on the visual centroid of the image because the scintillator 
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sensitivity is not included in the presentation of the raw data, however it was taken into account for 
computing the mean energies. Our results show that manipulation of electron energy distribution 
using the deformable mirror is somewhat restricted. The final result after optimization does not 
reach the objective mask completely despite that the mean energies can be varied by up to 10%. 
This result is somewhat unsurprising as the scope for controlling the electron spectrum is mostly 
limited by the physics of the interaction - while changing the transverse intensity profile can make 
big differences to the shape of the plasma electric field structure, changing the maximum field 
amplitude of the wakefield (and therefore peak energy of the accelerated electrons) will be limited. 

Numerical simulations Although the details of the initial conditions required for optimal beams 
are difficult to determine and are therefore found using the genetic algorithm, we can at least 
demonstrate how modifications to the phase front of the laser pulse can improve the beam proper¬ 
ties with an example. To illustrate the underlying physics of the plasma wave dynamics determined 
by the conditions of the driving laser pulse, we performed two dimensional (2D) particle-in-cell 
(PIC) simulations using the OSIRIS framework^^. Parameters similar to the experiment conditions 
were used, with a Gaussian plasma density profile to enable trapping of electrons in the density 
down ramp (see ref. 28 and Methods for details on the simulation). 

It was previously shown in ref. 29 that the focusing fields of laser plasma accelerators can be 
controlled by tailoring the transverse intensity profile of the laser pulse using higher-order modes, 
where generalization to 3D was also discussed. Here, we simulated a laser pulse with a fun¬ 
damental Gaussian mode (TEMqo) or a coherent superposition of a fundamental (TEMqo) and a 
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second-order Hermite-Gaussian (TEM 02 ) mode (Fig. 5a). Although the plasma wave has a larger 
amplitude when it is driven by a single mode laser pulse, the wake phase front evolves a backward 
curvature when electrons are trapped and accelerated (see top panel in Supplementary Movie 1 
and Fig. 5b). Contrastingly, the evolution of the wakefield driven by the laser pulse with additional 
mode forms a flatter plasma phase front at the point of trapping (Fig. 5d). 

In Fig. 5c and e, the momentum distribution of the forward accelerated electrons has a larger 
transverse spread for the single mode laser pulse compared to the one with the addition of higher 
order modes. This is a consequence of the different trapping conditions and accelerating fields 
from the coherent plasma wakefield structure, which is governed by the structure of the driving 
laser pulse. In a comparative test to show this effect is not simply due to a lower intensity, we 
repeated the simulation using a single fundamental Gaussian mode laser pulse with a larger focal 
spot with the same peak intensity as that with the superimposed modes. The wakefield evolution 
shows very similar response as Fig. 5b and does not develop a flatter phase front as seen in Fig. 5d. 
The subsequently accelerated electrons have very similar divergence to that in Fig. 5c, eliminating 
the possibility that the improvement comes from either a high intensity effect or a simple change 
in /-number. Note that we are not saying that this mixed-Hermite-Gaussian mode is the optimal 
pulse in the experiment; this is simply an illustration of how small changes in pulse shape can have 
significant effects on electron beam properties. 
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3 Discussion 


When a particular wavefront of laser light interacts with plasma, it can affect the plasma wave struc¬ 
tures and trapping conditions of the electrons in a complex way. For example, Raman forward scat¬ 
tering, envelope self-modulation, relativistic self-focusing, and relativistic self-phase modulation^*' 
and many other nonlinear interactions modify both the pulse envelope and phase as the pulse prop¬ 
agates, in a way that cannot be easily predicted and that subsequently dictates the formation of 
plasma waves. Moreover, under realistic experimental conditions, ionization dynamics before the 
laser pulse reaches the vacuum focus can also modify the phase of the driving pulse. Ideally, the 
light interacts in such a way as to generate large amplitude plasma waves with electric field struc¬ 
tures that accelerate electrons with small divergence, high charge etc. Because of the complicated 
interaction, it is difficult to determine a laser phase profile that will lead to such a plasma structure. 
However, such unforeseeable conditions were successfully revealed by using the evolutionary al¬ 
gorithm method, with the result that the electron charge can be increased and emitted in a very well 
collimated beam. 

Here we have implemented coherent control of a nonlinear plasma wave and demonstrated 
an order of magnitude improvement in the electron beam parameters. The laser beam optimized 
to generate the best electron beam was not the one with the ‘best’ focal spot. Control and shaping 
of the electron energy distribution was observed to be less effective, but was still possible. The 
capability for wavefront control was also limited by the number of actuators and maximum defor¬ 
mation of the deformable mirror used in our experiments. In addition, this work was performed 
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using adaptive optics, but it is clear that coherent control of plasma waves should be possible in a 
variety of configurations, for example by using an acousto-optic modulator to control the temporal 
phase of the driving pulse. Recently developed techniques^*’^^ for single-shot diagnosis of plasma 
wave structures may provide an avenue for direct control of the plasma evolution. 

The concept of coherent control for plasmas opens new possibilities for future laser-based 
accelerators. Although still at the stage of fundamental research, laser wakefield accelerators are 
showing significant promise. In principle, such improvements could be integrated into next gener¬ 
ation high-power laser projects, such as ICAN^^, based on coherent combination of many indepen¬ 
dent fibers, taking advantage of both their high repetition rate and controllability. The stability and 
response of the wakefield to laser conditions, such as phase front errors, is not well understood, 
but is crucial for the success of laser wakefield acceleration as a source of relativistic electrons and 
secondary radiation. For example, the presence of an asymmetric laser pulse was shown to affect 
the betatron oscillations and properties of x-rays produced in laser wakefield accelerators^"^"^’. Im¬ 
plementing the methods of this study should enable a significantly improved understanding and 
control of the wakefield acceleration process with regard to stability, dark current reduction and 
beam emittance. 

Methods 

Laser System The Relativistic Lambda Cubed laser (A^) produces 30 fs pulses of 800 nm light at 
a repetition rate of 500 Hz with an ASE (Amplified-Spontaneous-Emission) intensity contrast of 
~ 10® around 1 ns before the main pulse. The system is seeded by a Eemto-Laser Ti:sapphire os- 
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cillator, which generates 12 fs pulses and has a companion carrier envelope phase locking system. 
An RF addressable acousto-optic filter called a Dazzler controls the spectral amplitude and phase 
of these pulses. Selected pulses from the Dazzler train are stretched to 220 ps in a low-aberration 
stretcher and amplified to 7 mJ in a cryogenically cooled large-mode regenerative amplifier (Re¬ 
gen). The energy dumped from the Regen cavity is ‘cleaned’ in a Pockels cell and used to seed a 
3-pass amplifier as an upgrade from the laser system described in ref. 38, which delivers up to 28 
mJ pulses before compression. Following 71% efficient compression, 20 mJ pulses are trimmed to 
18 mJ at the perimeter of a 47 mm-diameter, 37-actuator deformable mirror. Throughout the sys¬ 
tem, pump light is provided by a variety of internally doubled Nd-doped YAG, YLF and vanadate 
lasers. The output beam with its controllable wavefront is then delivered to one of five experi¬ 
mental areas for the production of x-rays, electron beams, ion beams, THz radiation, high-order 
harmonics, or warm-dense matter. 

Electron Acceleration and Detection The focused laser pulse drives plasma waves by interacting 
with an argon gas jet flowing continuously from a 100 /rm inner diameter fused silica capillary. 
Typically the laser axis is 300 /rm above the orifice of the tubing. The laser pulses were focused 
by an //2 off-axis parabolic mirror to a spot size of 2.5 /um FWHM with a maximum of 10 mJ 
energy on target. The plasma electron density is measured to be in the range (0.5-2) x 10^® cm“^ 
using transverse interferometry. Electrons are accelerated in the density down ramp with a final 
energy in the 100 keV range^^ and detected by a high resolution scintillating screen (J6677 FOS 
by Hamamatsu), which is placed about 35 cm downstream from the source and imaged with a lens 
coupled 12-bit CCD camera for a 4x4 cm effective area. The scintillator sensitivity was calibrated 
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using an electron microscope for the energy range in the spectrum measurement. Electron beam 
charge was estimated using the calibrated scintillator response, manufacturer-provided information 
for the CCD camera (gain, quantum efficiency etc.) and the measured effective numerical aperture 
of the imaging system. 

Focal characterization and wavefront measurement The amplified laser beam was attenuated 
by using a half-wave plate and the polarization dependent properties of the compressor grating of 
the laser system. 025 mm neutral density filters (Thorlabs, Inc.) were inserted in the exit beam 
after the compressor and before a telescope beam expander. The laser focus was imaged by a 
60 X microscope objective lens (Newport Corporation, M-60X) onto a 8-bit CCD camera for focal 
characterization (cf setup A in Fig. 1). A Shack-Hartmann wavefront sensor (Flexible Optical 
BV) was used to determine the relative wavefront change between different deformable mirror 
configurations. The sensor, which consists of a 30 x 30 microlens array having a focal length 
of 3.5 mm and 150 /um pitch, was directly placed in the path of the converging beam after the 
focusing parabolic mirror. Reconstruction of the wavefront was performed using the FrontSurfer 
analysis software (Flexible Optical BV), typically with ~ 450 measured local wavefront slopes 
and RMS error on the order of 0.05A. Rotating the neutral density filters and the half-wave plate 
did not change the focal spot or the wavefront measurement significantly, insuring the wavefront 
distortion introduced by attenuation was negligible. 

Deformable Mirror and Genetic Algorithm The deformable mirror (AOA Xinetics) has a 47- 
mm clear aperture of a continuous face sheet with 37 piezoelectric actuators arranged on a square 
grid spaced 7 mm apart. The maximum stroke used in this experiment is about 2 fim. 
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The m irror shape is controlled by a genetic algorithm, which is a method mimicking the 
process of natural selection and routinely used to generate optimal solutions in complex systems 
with a large number of variables. The genetic representation in our experiments comprises a set 
of 37 independent voltage values for the deformable mirror actuator array. A fitness function is 
designed to produce a single figure of merit (FOM) to evaluate how close the solution is to the 
goal. 


In the electron beam profile optimization experiment, FOM is computed as follows: 
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where fj is the pixel intensity for every pixel {i,j) in the whole image and fq is a coordinate point 
in the image used as an optimization target. The power factor n > 0 gives higher weighting to 
those pixels closer to the target (inverse distance weighting). 


In the experiment to control the energy spectrum, FOM is calculated using the following 


formula given a pre-defined image mask. 


FOM = 1 - 
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The mean intensity is the sum of the pixel counts divided by the number of pixels for a defined 
region. A rectangular mask was used in the experiment as specified by the region enclosed by the 
red dashed lines in Fig. 4b. 


Numerical simulations The 2D PIC simulations were performed in a stationary box of the di¬ 
mensions 573 X 102 pm with 10000 x 600 cells and 4x4 particles-per-cell. A Gaussian plasma 


13 




density profile was used in the propagation dimension (xi), peaked at xi = 200 /rm with a full- 
width-at-half-maximum of 120 fxm and a maximum electron density of O.OOSric, where ric is the 
plasma critical density. The laser pulse was initialized at the left edge of the simulation window 
and focused at 215 /rm in the density down ramp. In 2D geometry, the transverse intensity pro¬ 
file of the laser pulse for fundamental Gaussian mode (TEMoo) has the form alexp{—2xywl), 
and the second-order Hermite-Gaussian mode exp(— 20 : 2 /( 8 ^ 2 “ 2)^, where ao ,2 is the 
normalized vector potential and wo ^2 is the beam waist parameter. The two modes are coherently 
superimposed in the same plane of polarization. Here we used even-order Hermite-Gaussian mode 
(TEM 02 ) for its symmetric property. A phase difference of tt/S was applied at the beam waist 
between the two modes to simulate variations in the optical phase front condition. The beam waist 
was positioned to account for focal shift as a result of coherent superposition of two modes such 
that the location of the maximum on-axis laser intensity was the same (at xi = 215 /um) for all 
simulation runs. The simulation parameters are ag = 1.08 and wg = 3.31 jum for the Gaussian 
mode alone, or ag = 1.0, 02 = 0.15 and wg = 102 = 3.31 jum for the superimposed mode. 
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Figure 1: Experimental set-up. DM: deformable mirror; OAP: off-axis parabolic mirror, 50 mm di¬ 
ameter; OL: objective lens; PD: photodiode. A: conventional focal spot optimization using second- 
harmonic generation; B: setup for direct optimization of the electron signal from the laser plasma 
accelerator. 
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Figure 2: Optimization of the eiectron spatiai profiie. Electron beam profile image integrated 


over 50 shots (100 ms exposure time) with a deformable mirror configuration (a) corrected for 


the best focal spot (BFS) and (b) 30 V on all actuators, (c)-(h) are single-shot electron beam 
profiles after genetic algorithm optimization using different weighting parameters, n. (i) shows the 


convergence of the genetic algorithm with n=8. The shaded gray area represents the range of the 


10 best children in each iteration and the solid green curve is the average, (j) Comparison of the 


peak charge density in a single-shot electron image. Error bars represent the root mean square 


shot-to-shot fluctuations; Estimated charge per shot is displayed for every image. Contours shown 


are for 20, 40, 60 mrad, centered on the beam centroid. 
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Figure 3; Laser intensity profiie in the focai region and wavefront. Scan of laser focal Intensity 
in vacuum with the deformable mirror optimized for (a) second harmonic signal (highest intensity) 
and (b) electron beam in Fig. 2. (c) Relative wavefront change reconstructed from direct mea¬ 
surement using a Shack-Hartmann wavefront sensor. Root-mean-square phase deviation in the 
aperture is 0.14 wave. The wavefront was reconstructed over a slighly smaller aperture (2.26 mm 
diameter) than the full beam diameter on the sensor (Ri2.7 mm 1/e^ width) to reduce errors in the 
peripheral area. (Scale bar, 1 mm). 
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Figure 4: Control of electron energy distributiion. (a) Schematic set-up for measuring electron 
energy distribution using dipole magnets, (b) Raw data showing the dispersed electron signal 
after genetic algorithm optimization using three different image masks. The location of the mask is 
indicated by the red rectangle and the black cross (x) represents the final mean energy for each 
spectrum. 50 shots were integrated for each spectrum. 


23 












jr/2 

jr/4 

VI 

0 « 
& 

-jr/4 


-jr/2 


electron density njnc 



Xi ijim) 

(jj njuc 



Xi Om) 


density 



PilmgC 


Figure 5; Modelling the effect of laser phase front using coherentiy superimposed modes. 

(a) Transverse laser intensity profile and phase at the focal plane (xi = 215 for a single 
Gaussian mode (solid curves) and a superimposed mode (dashed curves) obtained in the 2D 
simulations. Plotted are the laser pulses at vacuum focus. (b)(d) are snapshots showing the 
different plasma wave structures around the time {t = 1.28 ps) and location of trapping, driven by 
a laser pulse of single mode (b) and superimposed mode (d). The laser pulse propagates to the 
right. The phase space p 2 -pi distribution of the accelerated electrons shown for (c) single mode 
and (e) superimposed mode before the electrons exit the simulation box (t = 2.25 ps). 
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